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A Neptune Aerocapture Systems Analysis is completed to determine the feasibility, bene- 
fit and risk of an aeroshell aerocapture system for Neptune and to identify technology gaps 
and technology performance goals. The high fidelity systems analysis is completed by a five 
center NASA team and includes the following disciplines and analyses: science; mission de- 
sign; aeroshell configuration screening and definition; interplanetary navigation analyses; 
atmosphere modeling; computational fluid dynamics for aerodynamic performance and da- 
tabase definition; initial stability analyses; guidance development; atmospheric flight simula- 
tion; computational fluid dynamics and radiation analyses for aeroheating environment 
definition; thermal protection system design, concepts and sizing; mass properties; struc- 
tures; spacecraft design and packaging; and mass sensitivities. 

Results show that aerocapture can deliver 1.4 times more mass to Neptune orbit than an 
all-propulsive system for the same launch vehicle. In addition aerocapture results in a 3-4 
year reduction in trip time compared to all-propulsive systems. Aerocapture is feasible and 
performance is adequate for the Neptune aerocapture mission. Monte Carlo simulation re- 
sults show 100% successful capture for all cases including conservative assumptions on at- 
mosphere and navigation. Enabling technologies for this mission include TPS manufactur- 
ing; and aerothermodynamic methods and validation for determining coupled 3-D convec- 
tion, radiation and ablation aeroheating rates and loads, and the effects on surface recession. 


Symbols/Nomenclature 


A 

= Area (m 2 ) 

C^trim 

= Trim Angle of Attack 

CA 

= Axial Force Coefficient 

CBE 

= Current Best Estimate 

CD 

= Coefficient of Drag 

CFD 

= Computational Fluid Dynamics 

CG, eg 

= Center of Gravity 

CL 

= Coefficient of Lift 


CN 

= Normal Force Coefficient 

D 

= Drag 

GA 

= Gravity Assist 

L 

= Lift 

L/D 

= Lift-to-Drag ratio 

M/CDA 

= Ballistic Coefficient (kg/m 2 ) 

SEP 

= Solar Electric Propulsion 

TPS 

= Thermal Protection System 


I. Introduction 

A EROCAPTURE significantly increases the mass that 
can be delivered in orbit at a destination with an atmos- 
phere compared to an all-propulsive vehicle at the same des- 
tination with the same launch vehicle. Aerocapture utilizes 
aerodynamic forces on a vehicle during a single pass through 
a destinations atmosphere to capture into orbit about that 
destination, instead of a large propulsive delta V maneuver. 
An aerocapture flight profile schematic showing the primary 
aerocapture event sequence is shown in Fig. 1. * 1 

Aerocapture at Neptune is characterized by high entry 
velocities (28-30 km/sec inertial) into a H 2 He atmosphere, 
and capture into a high energy science orbit enabling Tritan 
flybys. Table 1 provides a comparison of the Neptune aero- 
capture reference mission, described in this paper, to a repre- 
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Figure 1. Aerocapture trajectory schematic. 
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sentative Mars aerocapture mission, and a Titan aerocapture reference mission 1 . The high entry velocities at Neptune 
compared to Titan and Mars result in significantly more severe environments at Neptune, including both aeroheating 
and g’s. The high energy science orbit for Neptune compared to the reference Titan and Mars missions, requires a 
significantly greater vehicle lift to drag ratio to provide adequate corridor width at Neptune. 

Table 1 Neptune aerocapture parameters compared to those at Titan and Mars. 



Neptune 

Titan 

Mars 

Entry Velocity (km/sec) 

29 

6.5 

5.7 

Nom. Entry Flight Path Angle (deg) 

-12.818 

-36 

-14.2 

Apoapsis/Science Orbit (km) 

3986 x 430,000* 

1700 

1400 

Atmosphere Composition (% volume) 

80% H2, 19% He, 
1% CH4 

95% N2, 

5% CH4 (max) 

95.3% C02, 
2.7 %N2 

Atmos Scale Height at Aerocapture Alt (km) 

49 

40 

10.5 

Atmospheric Interface Altitude (km) 

1000 (above 1 bar) 

1000 

125 

Aerocapture Altitude (km) 

100-300 (above 1 bar) 

200-400 

40 

Aerocapture Exit/Escape Velocity 

.97 

.69 

.76 

L/D 

.8 

.25 

.25 

M/CDA (kg/m 2 ) 

895 

90 

148 

Theoretical Corridor (deg) 

2.27 

3.5 

~1.4 

Time from Atmos Entry to Atmos Exit (min) 

10 

42 

10 

Convective Stag Point Heat Rate (W/cm 2 ) 

8000 

46 (.91 m nose rad) 

30 (1.9 m nose radius) 

Radiative Stag Point Heat Rate (W/cm 2 ) 

4000-8000 

93-280 

Negligible 

Max g’s During Aerocapture (Earth g’s) 

22 

3.5 

2.5-3 


* For set up of Triton flyby resonance at 488,000 or 393,000 km apoapsis 


II. Approach 

A multi-center aerocapture systems analysis team, including NASA engineers and scientists from Ames Re- 
search Center (ARC), the Jet Propulsion Laboratory (JPL), Johnson Space Center (JSC), Langley Research Center 
(LaRC), and Marshall Space Flight Center (MSFC), led by Langley Research Center, was kicked off in October 
2002 and completed in October 2003. The effort was funded through the Code S In Space program. 

The mission objectives and initial spacecraft design for the reference concepts are based on JPL’s TeamX study 2 
of the Neptune Orbiter with probes mission. From this starting point, further science definition and initial analyses 
are completed to provide understanding of the vehicle requirements and selection of the reference concept and mis- 
sion. Higher fidelity analyses are completed on the reference concept including mission design; aeroshell configura- 
tion screening and definition; interplanetary navigation analyses for determination of approach navigation delivery 
dispersions; atmosphere modeling; computational fluid dynamics (CFD) for aerodynamic performance and database 
definition; initial stability analyses; guidance development; atmospheric flight simulation; CFD and radiation analy- 
ses for aeroheating environments; TPS design, concepts and sizing; mass properties; aeroshell and spacecraft struc- 
tural design and sizing; spacecraft design and packaging; and mass sensitivities. 

A. Science 

The Neptune mission includes a Neptune orbiter and two probes. The orbiter science mission includes two years 
in Neptune orbit. The science orbit is selected to enable Triton flybys. The Neptune Orbiter science instruments 
were selected to be representative, and include visible imager, IR imaging spectrometer, UV imaging spectrometer, 
thermal-IR imaging spectrometer, ion and neutral mass spectrometer, magnetometer, charged-particle detector, 
plasma wave spectrometer, microwave radiometer, USO (radio occultations) and two identical probes. 3 
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B. Mission Design and Reference Concept 
Selection 

Many alternate mission designs are considered, 3 4 
including launch on Delta IVH and Atlas 551; grav- 
ity assists utilizing various combinations of Venus, 
Earth and Jupiter; SEP at various power levels or 
chemical stages; and aerocapture versus chemical 
insertion. Launch dates after 2015 are considered, to 
provide time for technology development. An SEP, 
aerocapture system is baselined for the reference 
architecture. 

The reference mission selected is a compromise 
between trip time, net delivered mass, inertial entry 
velocity, theoretical corridor width and aeroheating. 
Fig. 2 and 3 illustrate the net delivered mass and 
entry velocity vs. flight time for a range of 
SEP/aerocapture concepts. In general, as flight time 
decreases the net delivered mass decreases and entry 
velocity increases. The Delta IVH VJGA trajectories 
are selected for the reference mission concept based 
on delivered mass capability. Trip times less than 10 
years are eliminated due to the rapid decrease in de- 
livered mass capability and rapid increase in entry 
velocity (and corresponding aeroheating) with 
shorter trip times. 

To further select a trip time, entry velocity, and 
required vehicle L/D, an initial trade in available 
corridor width as a function of vehicle L/D and en- 
try velocity is completed. Fig. 4 shows the theoreti- 
cal corridor width vs. entry velocity and L/D for 
aerocapture into a 350,000 km apoapsis orbit at 
Neptune. The results in Fig. 4 illustrate several 
points. Theoretical corridor increases with both L/D 
and entry velocity. An L/D = 0.8 vehicle at 28 
km/sec provides approximately the same theoretical 
corridor width as an L/D = 0.6 vehicle at slightly 
greater than 30 km/sec. In addition, to achieve rea- 
sonable theoretical corridor widths for aerocapture 
into the high energy elliptic orbit requires vehicles 
with significantly greater L/D than the high heritage 
blunt body configurations. (Ex., The theoretical cor- 
ridor width is only approximately 0.8 degrees for an 
L/D = 0.25 for a 70° sphere cone.) 

To provide an initial estimated theoretical corri- 
dor required for comparison to the available corridor 
width, a quick corridor margin analysis is completed 
for combinations of L/D = 0.6, 0.8 and 28 and 30 
km/sec entry velocity. The theoretical corridor must 
be adequate to accommodate dispersions, uncertain- 
ties, and variability in approach navigated states at 
atmospheric interface, aerodynamics, atmosphere, 
and guidance robustness. Initial estimates for navi- 
gated errors show ±.5° 3a errors in entry flight path 
angle at 28 km/sec and ±.6° 3a errors in entry flight 
path angle at 30 km/sec. Aerodynamic uncertainties 
can be conservatively estimated to result in ±.2 L/D, 
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Figure 2. Delivered mass vs. trip time for a range of 
SEP/Aerocapture mission concepts considered. 4 
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Figure 3. Entry velocity vs. trip time for a range of 
SEP/Aerocapture mission concepts considered. 4 
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Figure 4. Theoretical corridor width available as a 
function of vehicle L/D and entry velocity. 
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guidance is estimated to capture 95% of theoretical corridor, and corridor loss due to total mean variability of the 
atmosphere is estimated using lift up and lift down trajectories at the global extremes of mean density in the initial 
NeptuneGRAM atmosphere model. Note that this approach does not account for high frequency variability in at- 
mospheric density and the corresponding impact on vehicle performance and margin. 

Results of this estimate are shown 
in Fig. 5 for the 4 cases considered, 

L/D = 0.6 and 28 km/sec, L/D = 0.6 
and 30 km/sec, L/D = 0.8 and 28 
km/sec, L/D = 0.8 and 30 km/sec. The 
L/D = 0.6 and 28 km/sec case show 
the estimated corridor loss due to ap- 
proach navigated errors, aerodynamics 
uncertainty, and total mean atmosphere 
variability and uncertainty. If these 
losses are RSS’d a required theoretical 
corridor width can be estimated for use 
in comparison to the total available 
theoretical corridor width. For this 
case, the required and available theoretical corridor widths are approximately the same, well within the accuracy of 
the estimate, and again, high frequency atmosphere perturbations are not included. Because of the early phase of 
design, the objective for selection of the reference concept is to select a concept with margin greater than the RSS 
required corridor width. As a result, an L/D = 0.8 vehicle with a 29 km/sec entry velocity is selected for the refer- 
ence concept, with an L/D = 0.6 vehicle kept as an option. 

The reference concept is therefore described as follows. The mission launches February 17, 2017, on a Delta IV 
H. The launch vehicle fairing is 5 m in diameter with a 4.572 m static payload diameter. The total launch capability 
is 5964 kg, with a launch C3 of 18.44 km2/sec 2 . The SEP system is a 30 kW EOL, 6 engine SEP system, that oper- 
ates to 3 AU. A Venus, Jupiter gravity assist is utilized. The total trip time is 10.25 years, with Neptune arrival in 
2027. Two probes are released at E-4 months (1 week apart). The probes enter at E-4 hours and E-2 hours. Aerocap- 
ture inertial entry velocity at Neptune is 29 km/sec, atmospheric interface is 1000 km above 1 bar. The orbit is 157° 
retrograde, 430,000 km by 3986km. The science mission includes two years in Neptune orbit for a total 12.25 year 
Neptune Orbiter mission. 

C. Aeroshell Configuration 

Aeroshell configuration screening is completed to develop shapes with L/D = 0.8, and L/D = 0.6 as an option, 
while maximizing volumetric efficiency and minimizing M/CDA. Several aeroshell shape classes are considered 
including ellipsleds, flattened ellipsleds, high fineness sphere cones, biconics and bent biconics, as shown in Fig. 6. 
The vehicle mass and volume are fixed, the shapes within each 
vehicle configuration class are varied parametrically. Newtonian 
aerodynamics, varified with CFD, is utilized to screen the con- 
figurations for L/D and M/CDA over a range of angle of attack. 

Packaging efficiency is also screened through determination of 
vehicle volumetric efficiency. Based on these analyses, the flat- 
tened ellipsled is selected for the reference aeroshell configura- 
tion. Details of the configuration screening and vehicle selection 
are contained in Ref. 5. 

D. Design Cycles 

Two design cycles are completed for the Neptune Orbiter. The 
original objectives were to package 3 probes within the aeroshell 
of the orbiter. As a result the design cycle one vehicle is 5.5 m in 
length with M/CDA estimated at 273 kg/m 3 , and an aeroheating 
design trajectory developed at 400 kg/m2 to provide mass growth 
margin. With this design, the system mass margin on the Delta 
IVH SEP VJGA was estimated to be less than the desired 35%. In 
addition, further analysis shows that the science objectives can be 
met with the probes carried externally and released prior to aero- 
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Figure 6. Configuration classes considered. 
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Figure 7. Comparison of Design Cycle 1 and Design Cycle 2 
configuration scale. 
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Figure 8. Delivery flight path angle dispersions at atmospheric 
interface for various navigation scenarios. 


capture. Two external probes are there- 
fore included in the design cycle two 
concept. The aeroshell is photographi- 
cally scaled from a 5.5 m length vehicle 
to a 2.88 m length, shown in Fig. 7. The 
entry allocation is 2200 kg, resulting in a 
ballistic coefficient of 895 kg/m 2 . This 
M/CDA is used for both the performance 
analyses and the design trajectories for 
the aeroheating and structure design. 


E. Navigation 

Navigation analyses are completed at 
28 and 30 km/sec with the Mars Recon- 
naissance Orbiter camera and an Entry - 
3 day data cut-off. In addition, results are 
completed with an MRO camera with 
two times the pictures and with an ad- 
vanced MRO camera. As a comparison 
results are also completed for an Entry - 
2 day data cut-off. Results are shown in 
Fig. 8. See Ref. 6 for detailed discussion 
on the navigation analysis and results. 

Monte Carlo simulations are com- 
pleted for the reference vehicle using 
+.51° 3a entry flight path angle disper- 
sions. This is approximately equivalent 
to an MRO camera with 2x pictures or to 
an advanced MRO camera, each with a 
more conservative Entry -3 day data cut- 
off. A change to Entry -2 day cut-off 
significantly reduces the delivery entry 
flight path angle dispersions as shown in Fig. 8. 

F. Atmosphere Modeling 

A NeptuneGRAM 7 atmosphere model is developed 
based on Voyager and other data. Variability includes all 
measurement uncertainty, residual uncertainty due to tur- 
bulence and waves, and the expected variability due to 
latitude, altitude, seasonal and time of day variations. The 
atmosphere composition is 80% H2, -19% He, -1% CH4. 
Fig. 9 illustrates the total mean density variability as a 
function of altitude. Note that the range of aerocapture 
altitudes is between approximately 100 and 300 km. The 
parameter Fminmax is utilized to define the range of den- 
sity profiles. The mean density profile is represented by 
Fminmax = 0, the minimum density profile is represented 
by Fminmax = -1, and the maximum density profile is 
represented by Fminmax = +1 . 

Fig. 10 illustrates the latitudinal variation of density for 
the particular arrival season of the reference concept. This 
variation of Fminmax with latitude is represented by 

Fminmax = 0.44*cos(4.0*latitude) + fbias 



Figure 9. NeptuneGRAM mean density profile 
variability. 8 


where -0.56 > fbias < 0.56 
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Fminmax 

Figure 10. Effect of latitude on NeptuneGRAM mean 
density profile for Neptune Orbiter arrival date. 8 
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Figure 11. High frequency density perturbations in 
NeptuneGRAM. 8 



Therefore, for a typical Neptune aerocapture trajectory flying through the atmosphere between 7-21° latitude, - 
0.6 > Fminmax < 0.93, compared to a range of -1.0 > Fminmax <1.0 for a global variation. 

Fig. 11 illustrates a sample high frequency density perturbation compared to the mean density. The mean density 
corresponds to a given Fminmax value in Fig. 9. rpscale controls the high frequency variability of the atmosphere, 
with rpscale = 1 representing the greatest expected variability for Neptune. The rpscale = 1 results in Fig. 9, show a 
sample of how the high frequency content can alter the mean variability. Note that the high frequency content can 
act to increase or decrease the mean density with altitude, in addition to adding the high frequency content. Rpscale 
= 0.5 represents a potential decreased high frequency content for the Neptune atmosphere. 

The reference concept performance, as shown below, is based on the latitudinal variation of Fminmax and 
rpscale = 1. Aeroheating and structure design trajectories are based on the full range of Fminmax and rpscale = 1. 


G. Aerodynamics 

The aerodynamic database is developed from viscous LAURA CFD analysis of the reference configuration. 5 The 
vehicle trims at 40° angle of attack with an axial eg location relative to the vehicle length of 0.51 aft of the nose, and 
a vertical eg relative to the vehicle length of .0166 below the vehicle waterline. For the trimmed vehicle 
L/D = 0.806, CD = 1.405, CL = 1.133. Initial stability analysis shows that the flat-bottom ellipsled is longitudinally 
and laterally stable. 5 

The aerodynamic uncertainties are based on the JSC el- 
lipsled analysis for Mars, consistent with the X-33 aerody- 
namic database uncertainty model in Ref. 9. As shown in 
Fig. 12, CA: ±0.048 and CN: ±0.12, each using base area 
as the reference. The trim angle of attack uncertainty is 
assumed to be ±4°, defined in this initial analysis to be 
double that for a typical blunt body, such as a 70° sphere 
cone. Cg uncertainties are ±0.5% for axial eg relative to 
the vehicle length and ±0.125% for radial cg relative to the 
vehicle length. Based on stacked aerodynamic uncertain- 
ties, the L/D uncertainty is ±26.4% and -22%. Based on an 
RSS of the aerodynamic uncertainties, the L/D uncertainty 
is ±13.5% and -14.3%. The Monte Carlo variability for 
2001 cases is between the RSS and stacked uncertainties. 9 

Note that the effects of large TPS recession and resul- 
tant shape change on the vehicle aerodynamics and cg lo- 
cation have not been quantified. This analysis was outside 
the study scope. 
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Figure 12. Aerodynamic uncertainties. 
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H. Guidance 

HYPAS guidance (ref. 10) was chosen for the Titan aerocapture systems analysis. HYPAS utilizes vehicle lift 
and bank angle control through the atmosphere to target the desired exit orbit apoapsis and inclination. It is an 
analytically derived algorithm based on deceleration due to drag and altitude rate error feedback. This analytic, non- 
iterative, on-the-fly approach leads to efficient code (-320 source lines in Fortran), minimal storage requirements, 
and fast and consistent execution times. 

HYPAS consists of two phases: 1) Capture Phase: Establishes pseudo-equilibrium glide conditions; 2) Exit 
Phase: Exit conditions are predicted analytically assuming a constant altitude rate followed by constant acceleration. 
The lift vector is adjusted to null the error between predicted and target apoapsis, and bank reversals are used to 
keep inclination errors within the desired limits. Results show excellent performance and an ability to capture -93% 
of the theoretical corridor. 


I. Performance/Simulation 

The reference concept performance is simulated in a Monte Carlo simulation 7 8 and includes each of the uncertain- 
ties and dispersions as described above. Fig. 13a-d show the reference concept Monte Carlo results. The reference 
concept is an L/D = 0.8 vehicle, with M/CD A = 895 kg/m 2 . The target orbit, to enable Triton flybys, is retrograde 
with an apaopsis of 430,000 km, and a periapsis of 3986 km. Uncertainties included in the Monte Carlo include 
navigation, with ±0.51° 3a entry flight path angle at atmospheric interface, atmosphere variability as a function of 
latitude and high frequency perturbations corresponding to rpscale = 1, and aerodynamic uncertainties described 
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Figure 13a. Reference concept Monte Carlo results, 
apoapsis vs. periapsis. 8 


Figure 13c. Reference concept Monte Carlo 
results, heat load vs. peak heat rate. 8 
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above. 

All 2001 Monte Carlo trajectories successfully capture. Approximately 66% of the cases achieve apoapsis within 
the apoapsis target range of 355,000 km and 500,000 km. In each of the Monte Carlo cases an apoapsis correction 
delta V, along with the periapsis raise delta V, is utilized to correct the orbit to the target of 430,000 by 3986 km. 
The total delta V, as shown in Fig. 13b, is 141 m/sec for the mean for the combined periapsis raise and apoapsis 
correction and 360 m/sec 99.87 percentile. Figures c and d illustrate the entry g loading and the peak heat rate vs. 
total heat load based on a stagnation point convective indicator for a lm nose radius. The 3a high g’s are 17.6 g’s, 
which are less than the 22.1 g’s used to design the vehicle structure. The 3a high heat rate and heat load stagnation 
point convective indicators are 2957 W/cm 2 and 295 kJ/cm 2 , respectively, compared to the 3250 W/cm 2 and 290 
kJ/cm 2 stagnation point convective rate and load indicators of 
the reference aeroheating design trajectory. 

The apoapsis error (prior to delta V correction) for Nep- 
tune is greater than that seen in previous studies. Before any 
apoapsis correction, and a 430,000 km apoapsis target, the 
3a range in Neptune apoapsis is 371,300 to 832,700 km. For 
comparison the range in apoapsis at Titan, 1 prior to any delta 
V to adjust apoapsis, and a 1700 km apoapsis target, is 1499 
km to 1883 km. The larger apoapsis errors at Neptune com- 
pared to Titan result from the high energy Neptune target 
orbit. At Neptune, the aerocapture exit velocity is very close 
to the Neptune escape velocity, resulting in a high sensitivity 
of apoapsis to aerocapture exit velocity. For example the 
ratio of the aerocapture exit velocity to escape velocity at 
Neptune is 0.97. The aerocapture to escape velocity at Titan 
is 0.69. 10 

The updated range of dispersions and uncertainties in 
navigation, aerodynamics, and atmosphere, are utilized in 
Fig. 13e, to assess the corridor margin for comparison to the 
original estimates in Fig. 5. The stacked aerodynamic uncer- 
tainties are used, and are similar to the assumptions earlier. Note that the atmosphere uncertainties are significantly 
less than the initial estimates. This results from the reduced range of Fminmax, by incorporating the variation of 
density with latitude as opposed to utilizing a global range, and also due to the higher vehicle ballistic coefficient 
and reduced atmosphere variability at lower altitudes. The revised estimates show significant margin above the RSS 
value. The effects of high frequency density perturbations and additional aerodynamic uncertainties due to surface 
recession are not represented in the estimate, however. Results suggest that there may be margin in the performance 
design that could be utilized to reduce the vehicle L/D requirement, reduce the entry velocity or to accommodate 
increased atmosphere variability resulting from an increase vehicle size (and lower M/CD A) in an effort to reduce 
aeroheating. 


J. Angle of Attack Modulation Option 

Utilizing angle of attack control as an option to augment 
the bank angle modulation is considered to assess any poten- 
tial benefits to performance and robustness. 10,8 Angle of at- 
tack modulation provides increased responsiveness to high 
frequency density perturbations and may assist with uncer- 
tainties in trim angle of attack. Angle of attack control could 
be provided with movement of an internal ballast or possibly 
with an aerodynamic control surface. Fig. 14a, b and c show 
results from the same Monte Carlo, one case without angle of 
attack control, and one case with ±5° angle of attack modula- 
tion. As shown, alpha modulation results in a significant re- 
duction in apoapsis dispersions, delta V and g’s. 

K. Aeroheating Environments 

The aeroheating design trajectory utilized for TPS sizing 
was based on the highest heat load trajectory from an earlier 
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Figure 14a. Monte Carlo results without angle of 
attack modulation. 8 
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Figure 14b. Monte Carlo results with angle of 
attack modulation. 8 
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Figure 14c. Comparison of performance parameters 
with and without angle of attack modulation. 


version of the reference Monte Carlo and is based 
on navigated uncertainties of ±.51° 3a, the global 
range of Fminmax variability, rpscale = 1, and the 
aerodynamic uncertainties described earlier. Fig. 
15a illustrates the range of heat rate vs. heat load for 
the Monte Carlo compared to the lift up lift down 
range. The Monte Carlo heat rate range is 2050- 
3250 W/cm 2 , and heat load range is 195-290 kJ/cm 2 . 
The lift up, lift down peak heat range is 3155-1122 
W/cm 2 , respectively. The lift up lift down heat load 
range is 185-442 kJ/cm 2 , respectively. Typically the 
vehicle is designed to fly significantly closer to the 
center of the lift up lift down heat rate and load 
range than shown for these results. In this case, the 
guidance is designed to fly lift down early in the 
entry trajectory to allow successful targeting of the 
high-energy orbit apoapsis with the high ballistic 
coefficient vehicle. Fig. 15b illustrates the time 
variation of the stagnation point heating indicator 
for trajectory #1647 compared to that for the mini- 
mum atmosphere lift up and maximum atmosphere 
lift down trajectories. This further illustrates that the 
design and corresponding Monte Carlo results are 
skewed toward the lift up high heat rate profiles. 

Because of these results, the peak heat load tra- 
jectory from the Monte Carlo, #1647, which also 
has -98 percentile peak heat rate of 2001 trajecto- 
ries, is selected as the reference trajectory for the 
TPS design, instead of the more traditional selection 
of the lift up trajectory for TPS selection, and lift 
down trajectory for TPS sizing. 

Turbulent convective (LAURA and DPLR) and 
radiative (NEQAIR and RADEQUIL) computations 
are completed on the reference vehicle 
(m/CdA= 895 kg/m 2 2.88 m flattened ellipsled) lift 
up and lift down trajectories and are utilized to es- 
timate “low”, “med”, and “high” aeroheating envi- 
ronments along Monte Carlo trajectory #1647. 11 
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Figure 15a. Aeroheating design trajectory Monte Carlo 
results for convective stagnation point heat load vs. heat 
rate on a 1 m nose radius. Illustration of aeroheating 
design trajectory #1647. 



Time (s) 

Figure 15b. Comparison of aeroheating profile for Monte 
Carlo trajectory #1647 to lift upminimum atmosphere and 
lift down maximum atmosphere trajectories. 
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Transition to turbulence prior to peak heating is expected due to significant ablation. Although only small differ- 
ences, less than 10%, result in comparisons of LAURA and DPLR laminar aeroheating rates, large differences occur 
in turbulent heating comparisons between LAURA and DPLR. The turbulence models used in the analyses included 
Cebeci- Smith algebraic or Wilcox k-Q model with LAURA; Baldwin-Lomax algebraic turbulence with DPLR. 
None of the turbulence models were developed for, or validated in, high Mach H2-He flows. 

Radiation is a significant contributor to the Neptune aeroheating environments. Both NEQAIR and RADEQUIL 
are utilized to estimate the radiative aeroheating environments. Significant differences between the two predictions 
result. To assist in understanding the aeroheating environmnents, analyses of Galileo are completed using NEQAIR 
and RADEQUIL for comparison with historical analysis and flight data, and for comparison to the Neptune Orbiter 
study vehicle. Current uncoupled analyses predict the same order of magnitude results, (between 45.4 kW/cm2 and 
78.5 kW/cm2 for the total uncoupled convection and radiation aeroheating) as the historical uncoupled analysis 
(63.3 kW/cm2 shown in Table 2). Engineering approximations, Galileo analysis and flight data indicate that the ef- 
fects of convection/radiation/ablation coupling must be 
considered. No tools exist for modeling convec- 
tion/radiation/ablation for coupled 3-D flowfields. 

(Galileo was modeled with 1-D assumptions.) Higher 
fidelity coupled models are expected to reduce the envi- 
ronments compared to uncoupled results. Development 
and validation of methods for modeling coupled con- 
vection/radiation/ablation 3-D flowfields is one of the 
technologies identified as enabling as a result of this 
study. 

Fig. 16 illustrates the division of the vehicle into 
zones, defined based on the vehicle structure, and to 
allow individual selection and sizing of TPS, based on 
the point with the highest rates and loads in each zone, to reduce overall TPS mass. Zone 1 and 2 comprise the heat- 
shield or forebody of the vehicle. Zones 3 and 4 including the base, comprise the backshell. The vehicle maximum 
diameter, also referred to as the waterline, occurs at the boundary of zone 2 and 4. 

Table 2. Comparison of Neptune Orbiter reference concept to Galileo. 



Figure 16. TPS zones. 



Galileo (Jupiter Dec 1995) 

Neptune Orbiter (study) 

Atmosphere composition 

86.2% H2, 13.6% He 

81% H2, 19% He 

Inertial entry velocity (km/sec) 

60 

29 

Atmos relative velocity (km/sec) 

48 

31.4 

Inertial FPA (deg) 

-6.835 

-12.818 

Trajectory 

Ballistic 

Lifting, guided, controlled 

Configuration 

44.25 deg sphere cone 

Flattened ellipsled 

Scale 

1.25 m diam (.291 m nose rad) 

2.88 m length 

M/CDA (kg/m 2 ) 

224, 229 

895 

Heat pulse duration 

~20 sec 

~200 sec 

Uncoupled stag pt peak heat rate 
(convec + radiative) (kW/cm 2 ) 

63.3 

16 

Coupled conv/rad/ablation (kW/cm 2 ) 

17.0 flight 12 , 28.0 analysis 13 

?? 

TPS stagnation point thickness (cm) 

14.6 

12.9 

TPS stagnation point recession (cm) 

4.6 

9.6 

TPS material - heatshield 

Nose piece: fabricated from billet of 
chopped molded carbon phenolic; 
tape-wrapped carbon phenolic flank 

Nose: carbon phenolic (manufacturing 
approach??); Wind: reduced density 
carbon phenolic (dev/ testing?) 


Figures 17a and b show the range “low”, “med”, and “high” of peak heat rate and load estimated based on the 
CFD and radiative aeroheating analysis for the highest heat rate location on both the vehicle nose and the vehicle 
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Figure 17a. Low, medium and high aeroheating 
results for zone 1. 



Fully Genesis Fully Genesis Reduced 

Dense Concept Dense Concept Density 

CP CP CP 


Figure 18a. TPS sizing results for zone 1. 14 



Fully Genesis Reduced Fully Genesis Reduced 

Dense Concept Density Dense Concept Density 

CP CP CP CP 


Figure 18b. TPS sizing results for zone 2. 14 


Figure 17b. Low, medium and high aeroheating 
results for zone 2. 

wind side for trajectory #1647. Note that after further 
analyses, including analysis of Galileo, the “high” esti- 
mate is well outside the expected range of aeroheating 
environments. These aeroheating environments are util- 
ized to complete TPS selection and sizing. 

TPS materials considered for the heatshield include 
carbon phenolic, reduced density carbon phenolic, and 
the Genesis carbon fiber form with carbon-carbon face 
sheet concept. Results of TPS sizing for the “Low” and 
“Medium” aeroheating are completed and shown in Fig. 
18a and b, for the nose and wind side, respectively. 14 
“Medium” levels are utilized for the Reference. The nose 
region is characterized by significant recession. Fabrica- 
tion of the tape-wrapped carbon phenolic or Genesis con- 
cept may not be possible for these environments. As a 
result a fully dense carbon phenolic is selected for the 
nose region of the reference concept. However, TPS 
thickness in the nose region is beyond current TPS manu- 
facturing experience for this shape and acreage. If the 
aeroheating rates and loads remain at the levels esti- 
mated, TPS manufacturing approaches will be enabling 
for the Neptune aerocapture mission. For the wind side, 
the reduced density carbon phenolic is selected, but addi- 
tional work is needed to design and assess the ability of 
this type of concept to accommodate the heat rates esti- 
mated. 


Zone 

Material 

Mass (kg) 

Zone 1 (Nose) 

Fully Dense CP 

204 

Zone 2 (Wind) 

Reduced Density CP 

293 

Zone 3 (Lee, Nose) 

PICA 

0.6 

Zone 4 (Lee, Nose) 

SLA 561 

58 


Figure 18c. Reference concept TPS selected and cor- 
responding CBE mass. 14 

L. Aeroshell Structure 

Fig. 19a and b show the reference vehicle structural 
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Backshell & Base 


Upper Thruster 
Support Frame 


Note: Deck will be 
trimmed in front of 
thrusters to facilital 
aeroshell separatic 

Deck Assembly 


Deck Support Ribs 

Integral Ribs 
and Longerons 


Foreboc 


Component 

CBE Mass (kg) 

Aeroshell 


94.85 

Forebody 

44.9 


Backshell 

42.8 


Base 

7.15 


Deck 


21.6 

Deck Ribs 


17.9 

Tank Supt Rods 


.42 

Thruster Supt Frame 


1.75 

Total Ellipsled Structure (CBE) 


136.5 


Figure 19a. Reference structural concept. 15 


Figure 19b. Reference concept CBE masses. 15 


concept and initial current best estimate of mass properties, respectively. 15 Optimization of the structure after com- 
pletion of the study indicates an opportunity to reduce the structural mass. 15 Launch loads and stiffness require- 
ments, and aerocapture entry loads are considered in the design and sizing of the structure. The load path for the 
orbiter on launch is from the cruise stage through the aeroshell to the deck. The TPS mass is considered to be a para- 
sitic mass. The aeroshell forebody, backshell, base and deck are 2.54 cm thick sandwhich construction with 5052 A1 
honecomb core and Graphite/Polymide face sheets. Integral monolithic blade stiffeners, longerons and ribs, are in- 
cluded for the forebody and backshell. The deck includes deck support ribs. 20 separation fittings attach the 
aeroshell forebody and backshell, and deck, which are used to separate the backshell and forebody from the deck 
after aerocapture. 


M. Mass Properties, Packaging 

Fig. 20a, b illustrate the packaging of the aerocapture orbiter, two probes and SEP propulsion module in the 5m 
Delta IVH fairing 3 . Fig. 20c illustrates the packaging of the aerocapture orbiter. Table 3 includes the mass summary 
of the reference vehicle concept. The stack wet launch allocation is 5500kg. The aerocapture entry allocation is 
2238kg (~2% greater than the allocation used in the performance analyses). 35% margin (allocation - 
CBE)/allocation is included on dry mass, with -8% unallocated launch reserve. The aerocapture mass fraction is 
59% of the orbiter dry mass based on growth masses (“MEV” in Table 3) with aerocapture propellant included 
(aeromaneuvering, periapsis raise and apoapsis correction); and 50% without aerocapture propellant included. 



Figure 20a. Reference 
concept packaging in Delta 
IV, 5 m fairing. 3 


Figure 20b. Detail of reference Figure 20c. Reference concept orbiter 
concept packaging in Delta IV, 5 m packaging. 3 
fairing. 3 
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Table 3. Reference concept mass property summary. 3 


Mass in kg 

UBE UonT MEV iviarg 


Launcn capability 


S9b4 

Launcn Keserve 

— 57T7o 

4b3 

Launcn wet aiioc 


55UU 

sth’ lv Adapter 

43 

3U.U% 

bZ 

12Z7o 

/u 

' Xenon 

y C6 

lU.UVo 

“TTT7TT 

U.U7o 

IU7U 

' stM ury iviass 

"TOT 

zy.bVo 

"TOF 

2U.U7o 

TFoZ 

' uruise hyarazine 


“ TTT 


TTT 

" uruise hrooes 

"TO 3U.lT/o 

2U/ 

2U.U7o 

2TF 

a/u tntry aiioc 


ZZ6ti 

a/u Aerosneii/ 1 hs 

73b 30/0% 


zU.U/o 

1149 

' A/U AUS Prop 


zz 


ZZ 

a/u peri Kaise prop 


“TO 


TO 

urmt wet aiioc 


“szf 

uroit prop 




T2T 

urmt ury Mass “ 

5Z4 7773% 

bb/ 

ZU.47o 

bU4 


Ubt = uurrent best tstimate 
Uont = Contingency = (MtV-Ubt)/Ubt 
MbV = Maximum txpected Value 
Marg = Margin = (Alloc-MtV)/MtV 
Alloc = Allocation 


N. Comparison to All-Propulsive Mission 

Several alternate mission concepts are shown in Table 4 4 for comparison to the reference concept labeled “Op- 
tion B2”. Each option shows the mass that can be delivered to Neptune prior to insertion, labeled “Pre-NOI Net 
Delivered Mass.,” and the mass required to capture into Neptune orbit. For the chemical insertion the chemical pro- 
pellant and chemical dry mass are calculated based on the “Pre-NOI Net Delivered Mass”. For the aerocapture sys- 
tem, the “Aerocapture System” mass is based on the reference concept and is fixed at 1119 kg. The “Payload in 
Neptune Orbit” is defined based on the reference concept and is 792 kg. “System Margin” represents either a sur- 
plus or deficit in the capability of the system to deliver the 792 kg into orbit. The System Margin should be between 
15-20% for adequate margin. 

To determine the benefit of aerocapture compared to an all-propulsive system, the aerocapture system that deliv- 
ers the maximum mass to Neptune orbit (Delta IVH, EJGA, SEP, Aero) can be compared to the all-propulsive sys- 
tem that delivers the maximum mass to Neptune orbit (Delta IVH, EJGA, SEP, Chem), each for the same launch 
vehicle. For the all-propulsive option, a maximum of 1167 kg can be delivered into Neptune orbit (zero margin). 
For the aerocapture option, assuming a fixed aerocapture mass fraction of 59% (includes aerocapture deltaV), 1614 
kg can be delivered into Neptune orbit (zero margin). Therefore aerocapture results in approximately 1.4 times 
more mass in Neptune orbit as compared to an all-propulsive system. 

In addition, Table 4 shows significant trip time savings for the aerocapture systems as compared to the all- 
propulsive systems. 
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Table 4. Comparison to alternate mission concepts. 4 


Launch Vehicle 

Delta IV H 

Atlas 551 

Gravity Assist 

VEJGA 

EJGA 

VJGA 

EJGA 

Earth to Neptune Prop System 

Chem 

Chem 

SEP 

SEP 

Chem 

SEP 

NOI Prop System 

Chem 

Aero 

Aero 

Chem 

Aero 

Chem 

Aero 

Aero 

Aero 

Option 

A1 

A2 

A2 

B1 

B2 

B1 

B2 

A2 

B2 

r 7 

Cruise Time to Neptune (yrs) 

15.0 

10.8 

11.8 

15.0 

10.5 

15.0 

10.3 

11.8 

10.5 

Launch Year 

2014 

2016 

2014 

2016 

2016 

2017 

2017 

2014 

2016 

Launch C3 (km2/sec2) 

15.6 

26.0 

47.3 

13.5 

13.6 

17.0 

18.4 

47.3 

9.1 

SEP Power (kW, EOL) 




30 

30 

30 

30 


30 

Inertial Entry Velocity (km/s) 


29 

29 


29 


29 

29 

29 

Neptune Cruise Chem DV (m/s) 1 

3429 

1413 

357 





357 


NOI Chem DV (m/s) 1 

2300 



2871 


2781 





r 9 

r 6 

6 6 

Launch Capability 

7012 

5695 

3550 

6543 

6532 

6130 

5964 

2630 

4850 

Propellant Mass?' 3 

4158 

2040 

376 

655 

809 

1025 

1070 

279 

713 

LV to Prop Module Adapter 

62 

62 

62 

62 

62 

62 

62 

62 

62 

Prop Module Dry Mass 

806 

542 

289 

1437 

1449 

1465 

1468 

243 

1441 

Chem Prop Mod to Payload Adapter 

40 

40 

40 





40 


Pre-NOI Separated Mass 10 

318 

318 

318 

318 

318 

318 

318 

318 

318 

Pre-NOI Net Delivered Mass 

1628 

2694 

2464 

4071 

3895 

3260 

3046 

1688 

2315 


Aerocapture System 4 


1119 

1119 


1119 


1119 

1119 

1119 

NOI Chem Propellant Mass 3 

966 



2417 


1898 




NOI Chem Dry Mass 

280 



487 


413 




Payload in Neptune Orbit 

792 

792 

792 

792 

792 

792 

792 

792 

792 

System Margin = LV-MEV 

(409) 

783 

553 

375 

1984 

157 

1135 

(223) 

404 

System Margin % = (LV-MEV)/MEV 

-5.5% 

15.9% 

18.5% 

6.1% 

43.6% 

2.6% 

23.5% 

-7.8% 

9.1% 


MEV: Maximum Expected Value = best estimate + 30% contingency 

Assumptions and Notes: 

All masses are MEV mass listed in kg 

1 Includes 5% DV contingency 

2 Chem Propellant mass calculated using "Launch Capability" as system total mass; Chem Isp = 325 sec 

3 SEP Propellant mass calculated using "Launch Capability" as system total mass; includes 10% prop mass contingency 
4 Aerocapture System Mass: aeroshell structure, TPS, and DVto achieve 28766x488,000 km orbit 

6 Propellant mass and Prop Module Dry Mass for SEP / Chem options includes propellant and dry mass for 

both SEP and chemical stages 

7 Neptune Aerocapture Study Reference Mission 

s Chem Propellant mass calculated using "Pre-NOI Net Delivered Mass" as Initial mass; Chem Isp = 325 

9 Total Cruise+NOI DV split equally between two stages; l.e. Cruise delta-V is staged 

10 Includes Probes and ~100kg of cruise hydrazine 


O. Summary and Technology 

Aerocapture can deliver 1 .4 times more mass to Neptune than an all-propulsive system for the same launch vehi- 
cle. Aerocapture is feasible and performance is adequate for the Neptune aerocapture mission. Monte Carlo simula- 
tion results show 100% success for all cases including conservative assumptions on atmosphere and navigation. Ad- 
ditional analyses are required to assess the amount of surface recession from coupled 3-D convec- 
tive/radiative/ablation analyses, determine the aerodynamics and uncertainties resulting from time and path depend- 
ent shape change, and evaluate the effect on guidance and control algorithm design, and performance. The Neptune 
spacecraft can be successfully packaged in an aeroshell and result in ~8% unallocated mass while meeting the re- 
quired mass margins. 

Technologies identified in the study as requiring development are grouped into three categories; enabling tech- 
nologies, strongly enhancing technologies and enhancing technologies. Technologies annotated with an asterisk are 
categorized based on current understanding. Additional assessment could change the categories. 

The enabling technologies identified include 

• TPS Manufacturing. TPS thicknesses are beyond current manufacturing experience for carbon phenolic for 
this shape and acreage. 

• Aerothermodynamic methods and validation 
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o Aerothermodynamics are characterized by high radiative and convective aeroheating, coupled convec- 
tion/radiation/ablation, and significant surface recession with effects on vehicle aerodynamics on a more 
complex shape. 

o Coupled convection/radiation/ablation capability for three-dimensional flowfields is needed for definition 
of aeroheating environments, TPS requirements, and vehicle shape change, 
o An approach is needed to determine and represent the aerodynamics/uncertainties on the time varying path 
dependent shapes and corresponding masses in an aerodatabase and simulation. 

The strongly enhancing technologies identified include 

• Guidance Algorithm* - Existing guidance algorithms have been demonstrated to provide adequate perform- 
ance. However, improvements are possible to improve performance, to determine the ability to reduce heat 
loads and to accommodate time varying, path dependent shape and ballistic coefficient change 

• Flight Control Algorithm* - Algorithms must be able to accommodate shape change uncertainties 

• Atmosphere Modeling - Neptune General Circulation Model output is needed to represent the dynamic vari- 
ability of the atmosphere. 

• Reduced Mass TPS concepts, ex., reduced density carbon phenolic, could be utilized to decrease aeroshell 
mass. 

• Utilizing the TPS as a structural element may reduce the combined structure plus TPS mass. 

• Alpha Modulation* reduces the dispersions in apoapsis, provides additional and more rapid response to density 
perturbations, and provides additional margin for trim angle of attack uncertainties. 

• Dual Stage MMRTGs 

• Deployable Ka-Band HGA 

The enhancing technologies identified include 

• Automated navigation, improved optical navigation camera. 

• Miniaturized ACS components. 

• Lower Mass, Power Science Instruments 


III. Future Work 

Several areas are recommended for future systems analysis in addition to the specific technology items listed 
above. Recommendations are as follows. 

Complete partial design cycles for one or more intermediate (between 2.88 m-5.5 m length scale) vehicle sizes. 
There may be a minimum mass vehicle between the 2.88 m and 5.5 m length vehicle. The trade is surface area vs. 
areal density of the combined TPS and structure. 

The current design has 460 kg unallocated mass. In addition, interplanetary trajectory designs have resulted in 
increased delivered mass capability. Several design changes can be considered within the increased mass capability. 
For example, an increased vehicle scale may reduce aeroheating rates and loads and the corresponding surface re- 
cession and TPS thickness required. 

Additional trades that can be completed include a further assessment of chemical vs. SEP cruise; additional sys- 
tems analysis considering angle of attack modulation; revisiting the L/D=.6 vehicle; consideration of a symmetric 
version of the flattened ellipsled i.e. an elliptic upper section, in addition to elliptic lower and section; utilization of 
the TPS as a structural element; and consideration of variable thickness TPS for TPS mass reduction. 
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